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ABSTRACT 

Polarization data from high mass star formation regions (W51 e2/e8, Orion 
BN/KL) are used to derive statistical properties of the plane of sky projected magnetic 
field. Structure function and auto-correlation function are calculated for observations 
with various resolutions from the BIMA and SMA interferometers, covering a range in 
physical scales from ~ 70 mpc to ~ 2. 1 mpc. Results for the magnetic field turbulent 
dispersion, its turbulent to mean field strength ratio and the large-scale polarization 
angle correlation length are presented as a function of the physical scale at the star 
formation sites. Power law scaling relations emerge for some of these physical quan- 
tities. The turbulent to mean field strength ratio is found to be close to constant over 
the sampled observing range, with a hint of a decrease toward smaller scales, indicat- 
ing that the role of magnetic field and turbulence is evolving with physical scale. A 
statistical method is proposed to separate large and small scale correlations from an 
initial ensemble of polarization segments. This also leads to a definition of a turbulent 
polarization angle correlation length. 

Subject headings: ISM: clouds — ISM: magnetic fields, polarization, turbulence — 
ISM: individual (W51 e2/e8, Orion BN/KL) — Methods: statistical 



1. Introduction 

Giant molecular clouds - the sites of star formation - are threaded by magnetic fields. The 
exact role of the magnetic field, the nature of turbulence and their interplay are still a matter of 
debate in the lite rature. Evidence for a weak magnetic field (super- Alfvenic turbulence) has been 



presented in e.g. ICrutcher et al.1 (I2009|) : IPadoan et al.l (|2004f) . whereas support in favor of a strong 
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magnetic field (sub-Alfve nic turbu l ence) controlling the formation and evolution of the molecular 
cloud is discussed in e.g. Ili et al.1 ([2009). Accurate measurements of the magnetic field strength 
are a key in distinguishing between these two theories. 

Different techniques have been employed to measure the magnetic field intensity and struc- 
ture on various scales. The Zeeman effect provides the only known method of directly measuring 
magnetic field strengths along the line of sight in a molecular cloud. Generally, it has to rely 
on strong enoug h line intensities and a l so high spectra l resolution in order to detect the splitting 



of sp ectral lines (IGoodman et all 1 19891 : ICrutcherlll999l : iTroland & Crutcheill2008l : ICrutcher et al. 



2009). Polarization of dust thermal emission at infrared and submillimeter wav elengths provides 
another method to study magnetic field properties (e.g.. iHildebrand et al.ll2000|) . The dust grains 
are thought to be aligned with their shorter axis parallel to the magne tic field lines, therefore the 
emitted light appears to be polarized perpendicular to the field lines (|Goldreich & Kylafislll98lL 



1982l : braine & Weingartneril 199611 1 9971: iLazarianl 2000 ). Radiative torques a re likely to be respon 



sible for the dust alignment ( Cho & Lazarian 2005; Lazarian & Hoan g 2007). Complementary to 



rqi 
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Zeeman splitting, dust polarization measurements are probing the plane of sky projected magnetic 
field direction. However, in order to derive the actual magnetic field strength perpendicular to the 
line of sight, additio nal assumptions are needed, as e.g. in the commonly used Chandrasekhar- 
Fermi (CF) method dchandrasekhar & Fermi 1953 ) or in its variations (e.g., Houde et al. 20041: 
Curran & Chrvsostomoull2007 ). 



When applying any of these CF methods, the dispersion in the measured plane of sky polar- 
ization is a key parameter. Most stud ies up to now relied on a magnetic field dispersion measurec 



2009a) or a model field (e.g., iGirart et al.l 120061 : iRao et al 



about a large scale mean field (e.g., IC hrysostom ou et al.1 1 19941: lLai et al.l 2001, 



2 009^ As noted in 



2002; Tang et al. 



Hildebrand et al. 



(120091) . in dense clouds the magnetic field structure might be the combined result from a va- 
riety of effects, such as differential rotation, gravitational collapse and expanding HII regions. 
Consequently, a globally derived dispersion might not reflect t he true contribution fro m magne- 
tohydrodynamic waves and/or turbulence. The recent work by IHildebrand et al.1 d2009|) develops 
a method based on a dispersion function about local mean magnetic fields. Besides providing a 
measure for the turbulent dispersion, the method also gives an accurate estimate of the turbulent 
to mean ma gnetic field streng t h ratio . Furthermore, the method is independent of any large scale 
field model. IHildebrand et al.l (|2009n discuss applications to the Orion, M17 and DR21 molecular 
clouds, observed with the Hertz polarimeter at the Caltech Submillimeter Observatory (CSO) with 
a resolution of 20". 

Dust polarization obs ervations have been carried out o ver a range of scales: fr om the large 



scale cloud envelope (e.g ., ISchleuningll 19981 : lLai et al.ll2001|) to collapsing cores (e.g jGirart et al 



20061 : iTang et al.ll2009bf) . The goal of this paper is to apply and extend the method developed in 
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Hildebrand et al.l (I2009Q across a range of scales in the same star formation regions. 



The paper is organized as follows: Section [2] gives a brief summary of the W51 e2/e8 and 
Orion BN/KL high mass star formation sites with the results relevant for our analysis. Section 
[3] defines the structure function and the auto-correlation function with some physical quantities 
resulting from the statistical analysis. The results are presented in section |4] with a discussion in 
section [51 Summary and conclusion are given in section [6] 



2. Data Set and Source Descriptions 



The presented data were obtained with the BIMA and SMA interferometers at wavelengths 
where the polarization of dust t hermal emission is traced. The detailed desc riptions of t he data 
analysis and images are given in Lai et al] ( 2001 ) for W51 e2/e8 with BIMA, Tang et ah d2009b ) 



for W51 e2/e8 with the SMA, iRao et al 



1998b for Orion BN/KL with BIMA and iTang et al. 



(|2010 ) for Orion BN/KL with the SMA. Relevant observation numbers are given in Tabled] The 
re-constructed features depend on the range of uv sampling and weighting. Nevertheless, these 
data currently provide the highest angular resolution (6>) information on the morphology of the 
magnetic field in the plane of sky obtained with the emitted polarized light in tho se star formation 



sites . They are thus complementary to the earlier analysis with single dish data in lHildebrand et al 



(|2009|) and references therein. This study is part of the program on the SMA0 (Ho, Moran and Lo 
2004) to study the structure of the magnetic field at high spatial resolutions. 



2.1. W51 e2/e8 



W51 e2/e8 are some of the strongest mm/s ubmm continuum sources in the W51 region. Lo- 
cated at a distance of 7 kpc (|Genzel et al.lll981|) . 1" is equivalent to ~ 0.03 pc. At a scale of 5 pc , 
measured at 850/xm w ith a resolution 9 « 15" with SCUBA on JCMT (|Chrysostomou et al.ll2002l : 
Matthews et al.ll2009|) . the polarization across the molecular cloud appears to be little organized 
and not uniform. However, the 350/im data, obtained with Her tz on the Caltech Submillimeter 
Observatory (CSO) sho w a well-ordered field (IDotson et al.ll2010h . A comparison of the Hertz and 
SCUBA data is given in Vail lancourt & Matthews! (I2008ab . Both the e2 and e8 core are unresolved 
by Hertz and SCUBA. lLaietal.1 (12001) reported a higher angular resolution (3") polarization map 
at 1.3 mm obtained with BIMA, which resolved out large-scale structures. In contrast to the larger 



1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and the Academia Sinica. 



-4- 



scale polarization map, the polarization appears to be more uniform across the envelope at a scale 
of 0.5 pc, enclosing the sources e2 and e8. In the highest angular resolution map obtained with the 
SMA at 870 txm with 9 O'.'J , the polarization patterns appear to be pinched in e2 and also pos- 
sibly in e8 (|Tang et al.ll2009b() . In the following, where not explicitly written in full, W51 always 
refers to W51 e2/e8. 



2.2. Orion BN/KL 



The Orion Mol ecular Cloud (OM C-1) is one of the closest massive star formation sites. At a 
distance of 480 pc (|Genzel et al.lll981|) . 1" is equal to 2.3 mpc. Source BN/KL is located near the 
strongest mm continuum emission, where the star formation pr ocess is active. Based on si ngle dish 
polarization measurements obtained with SHARP on the CSO. IVaillancourt et al.l (|2008b|) reported 
polarization maps observed at 350 and 450 fim in the OMC-1 ridge. The revealed magnetic fields 
appear to be uniform across the dust ridge with effective polarimetric beam sizes of 13" at both 
350 a nd 450 fim. Similar uniform polariz a tion m aps have been reported at 100 fim by ISchleuning 
dl998h and at 850 fjm bv Ivallee & Fieeel (12007]) with 9 « 15". Observed with BIMA at 1.3 mm 
and 3 mm with 9 of 3 "4 and 7". lRao et al.1 (11998) reported that the polarization appears to change 
abruptly in the south of the mm continuum peak. With the highest angular resolution achieved 
with the SMA at 870 /mi, 9 ~ 1", the polarization is consistent with the BI MA detections at 1.3 
mm and 3 mm, but the field appears to vary smoothly across the entire core (ITang et al.ll2010|) . In 
the following, where not explicitly written in full, Orion always refers to Orion BN/KL. 



3. Method 

This section summarizes how the statistical quantities are derived from the polarization data. 



3.1. Structure Function 



The polarization position angle (PA) structure function (of second order) is a measure of the 
mean square deviation in the plane of s ky projected magnetic field direction^ as a function of 
scale size. Following the recent work by lHildebrand et al.1 (120091) . we adopt their definition for the 
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magnetic field dispersion A(f>, which is the square root of the structure function: 



< A0 2 (4) > 



1/2, 



N(l k ) 



{Urd-Urjff 



1/2 



(1) 



h< r ij<h+\ 



where fa = PAj at position r t = (x^yt) and r ; y = a/ (Xj-Xj) 2 + (y,- -y ; -) 2 in each source reference frame 
with coordinates (xi,yi). < ... > denotes the averaging process over the entire polarization map 
with respect to r,y for each binning interval 4- N(lk) is the number of pairs of PAs with a separation 
in between the bins 4 and 4+i- All variables are plane of sky projected quantities. 



As derived in iHildebrand et al.l (I2009|) . assuming the magnetic field B to be composed of a 
smoothly varying large-scale mean field B and a statistically independent turbulent component B t , 
the structure function in the range 5 < l k <^d can be written as: 



< A0 2 (4) > tot ~ b 2 + m 2 ll + a 2 M (l k ), 



(2) 



where 6 is the turbulent field correlation length and d is the typical scale for variations in B . a Milk) 
are the binned error bars resulting from propagating the individual measurement uncertainties, b is 
interpreted as the scale-independent turbulent field dispersion and ml k is the linear dispersion term 
(with slope m) from the large-scale field B . All three contributions, being statistically independent, 
are added in quadrature. This is basically a first order Taylor expansion of the structure function 
in the domain where the turbulent field component, B t , is a s mall perturbation of a l arge-scale 
field B which is smooth on the scale of d. Further following IHildebrand et al.1 (|2009|) . the ratio 
between the turbulent and mean magnetic field strength can be calculated by evaluating explicitly 
the dispersion in the field directions under the assumptions of small perturbations, which results 
in: 

< B 2 t > x > 2 b 



V2-b 2 ' 



(3) 



3.2. Auto- Correlation Function 

The polarization angle correlation function measures the resemblance or self similarity (corre- 
lation with itself) of the projected polarization structure on average as a function of separation. Ad- 
ditionally, it leads to a definition of a characteristic polarization angle correlation length. Whereas 
the ratio in equation © relies on the assumption that the dispersion function can be written as a 
first order Taylor expansion on the smallest scales, the auto-correlation function with its weighted 
moments (e.g. correlation length in equation Q) is independent of such assumptions. In princi- 
ple, this provides an independent cross-check (see section I5TI) which can probe the assumptions 
in section [3TTI Furthermore, higher resolution data will allow for an even more detailed modeling 
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of the field structure without relying on a fit. We calculate the plane of sky projected polarization 
angle correlation function C as: 

<C(h)> = ^ £ Hrd-Hn) 

h<rij<h+\ 

= <0(r)-0(r+4)> (4) 

where the notation is identical to the one in equation ([I]). In order to make proper use of the auto- 
correlation function, one has to assume homogeneous isotropic turbulence, i.e. < C(4) > must 
depend only on the separations r (J . This requires 0, to be rotationally invariant. Whereas this is 
naturally the case for the structure function, which involves only the square of the difference of the 
PA, the product <^(r,-) • (j)j(rj) for the auto-correlation function depends on the reference frame and 
the definition of the range of the PA. The transformation — > £, (pj — > £— A0 ( y expresses all the 
correlation products with respect to the same position angle £(^0), where A&j = Since 
< A(pij < 90°, imposing zero correlation for perpendicular PAs fixes £ = 90°. This definition also 
ensures the correlation coefficients to be in the range between zero and one (PAs parallel) when 
normalized with £ 2 . 

Each observed PA will be the result of a superposition of a large-scale (</> ) and a turbulent 
contribution (5(f)). Consequently, the correlation function as written in equation dU) contains both 
contributions mixed. In order to characterize the large-scale polarization the turbulent part needs 
to be separated. Appendix A gives the details of how to derive a large-scale correlation function, 
< Co(h) >, from an ensemble of measured position angles </>, assuming <p = <po + 5(p. 

The characteristic large-scale polarization angle correlation length Ao is then calculated by 
integrating the weighted large-scale polarization angle correlation function: 

^ _ / < (f>o( r ) ■ M r +h) >r -h dh „ 
J < Mr) ■ <Po(r+h) >r dl k 
where the integration extends over the entire binning range. This is again a plane of sky projected 
quantity. In the case of a uniform polarization, with all PAs being parallel and aligned with a 
single direction, all the correlation coefficients will be one, independent of scale. The large-scale 
correlation length A is then in the middle of the largest and smallest scale, because all the scales 
are equally weighted (identical correlation) in this case. 

We note that the correlation length, being an integrated and weighted measure, is less sensitive 
to irregularities and incompleteness in the data set (see discussion in section 15711) . In analogy to 
equation ©, a turbulent polarization angle correlation length A, can be estimated, once large- 
scale and turbulent contributions are separated (Appendix A). A method to calculate the turbulent 
mag netic field correl ation length, based on a generalization of the dispersion function, is derived 



in 



Houde et all (I2009h . 
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4. Results 



Dispersion function and large-scale auto-correlation function together with turbulent corre- 
lation function are presented in Figure Q] and [2] for W5 1 e2/e8 a nd Orion BN/KL. For both high 
mass star forma tion regions, the pol arization data from BIMA (IRao et al] 1 19981 : lLai et al.ll200l|) 
and the SMA (|Tang et al.ll2009bll2010h are analyzed following section [3j Table Q] summarizes the 
observations and our findings. The binning intervals l k are set to integers of the synthesized beams. 
For elliptical beams, resulting from a non-uniform wv-coverage, the geometrical mean is adopted. 
Within each binning interval k, dispersion and auto-correlation are evaluated for 4 < r t j < 4+i • Cor- 
related data points below the synthesized beam resolution are removed. Only PAs with a polarized 
flux of more than 3<Ji p , the rms noise of the polarized intensity, are included. Errors of individual 
PAs are typically in the range of 5° to 10°. The binned error bars (om(4) in equation ©) in the 
Figures Q] and [2] are then determined by propagating the individual errors through the equations (OQ) 
and (HI). For the dispersion function they are typically around 0.°5 or less for the smallest scales. 
This is due to the sample variance factor from the large number of data points (~ 100 or more pairs 
of PAs). They grow to a few degrees at the largest scales. Fitting for the turbulent dispersion b is 
based on a least square fit of equation ([2]). The small binned error bars are neglected here, which 
is justified by the possible larger biases as it is discussed in section I5TT1 



4.1. W51 e2/e8 



The statistical analysis for W51 e2/e8 was performed with three data sets: BIMA (ILai et al. 



200 lb at 1.3 mm with a resolution of 2"3, covering the large-scale structure over ~ 20", and 
SMA at 870/xm w ith a resolution of 0"7, separately resolving the regions e2 and e8 over about 4" 
(|Tang et al. Il2009bh . All dispersion functions (structure functions of second order) show an increase 
over at least the first two bins (Figured! left panels). The larger scale BIMA measurement reveals 
a gentle increase in dispersion with a hint of a plateau at the larg est scales (Figure [H left bottom 
panel). This is very similar to the results in lHildebrand et alj (I2009|) . obtained with a 20" resolution 
in M17, DR21 Main and OMC-1. The higher resolution SMA observations show a steeper slope 
over the first two bins, with a dispersion at the smallest scales of about 40° and 55°, compared to 
about 10° in the BIMA observation. Whereas increase at smaller scales and tendency of a plateau 
at larger scales are still observed, the higher resolution observations show more irregularities. This 
is particularly the case for W51 e8. 



Following iHildebrand et al.1 (I2009|) the turbulent field dispersion b, as defined in equation ©, 
is obtained from the zero intercept of the fit at scale= 0. In order to stay in the linear regime, 
the first three bins from the BIMA data and only the first two bins from the SMA data are used. 



-8- 



(red solid lines in Figured] left panels). The resulting turbulent dispersions around the mean local 
magnetic field range from ~ 6° to ~ 54°, with corresponding turbulent to mean field strength ratios 
from 0.1 to 0.9 (equation ©). Higher resolution observations reveal larger values. 

The large-scale polarization angle correlation function for the BIMA observation shows a 
smooth curve as expected from the dispersion (Figure [Q bottom right panel, solid line): a small 
dispersion at small scales translates into a close correlation at these scales. At the BIMA largest 
scales, tracing the large-scale polarization variations, the auto-correlation decreases accordingly. 
The SMA observations of both e2 and e8 show correlations at the shortest scales which are fol- 
lowed by a rather sharp drop and a plateau-like extension. This again reflects the corresponding 
features in the dispersion functions. Both cores show a secondary peak in the auto-correlation 
function at larger scales, probably tracing symmetry features in the hourglass-like pinched field 
morphology. Calculating the characteristic polarization angle correlation length over the maxi- 
mum scale range, as introduced in equation ([5]), we find Ao = 230 mpc for the BIMA observation 
and Ao = 73 mpc and 63 mpc for e2 and e8, respectively. Due to the relatively small field of view 
sampled in our observations, the values of Ao possibly represent lower limits^ . This can then also 
explain why the correlation does not fully vanish at the largest scales. 

The same panels in Figure [TJ also show the small-scale correlations, separated by the method 
described in Appendix A. It is apparent that the turbulent correlation function, after an initial sharp 
drop, is still showing features similar to the large-scale function < C >■ This is a consequence of 
the weighting scheme outlined in Appendix A, where even at a larger scale small dispersion values 
can be accounted for a turbulent correlation with a certain probability. Limiting the turbulent 
correlation to within the first few bins yields A, between 25 and 45 mpc. 



4.2. Orion BN/KL 



Four data sets with very different resolutions were an alyzed for Orion BN/KL: BIMA observa- 
tions at 3 mm and 1.3 mm with a resolution of 7" and 3 ."4 ( Rao et allll998 ). and SMA observations 
at 870/mi from a combined compact with subcom pact, and a comp act with extended configuration 
with resulting synthesized beams of 2"8 and 0"9 (|Tang et al.ll2010|) . The general tendencies found 
for W5 1 e2/e8 - increase at smaller scales and plateau-like extension with some irregularities at 



Additionally, some information on the largest scales might be absent due to the missing zero-spacing in the 
interferometric observations. However, most of the observations used here contain information from short baselines 
with lengths comparable to a few antenna diameters, and coherent large-scale structures are apparent in the polarization 
maps. Correlations on even larger scales - also given the observed trend of decreasing correlation coefficients with 
larger scales - are therefore likely to add only negligibly to Aq. 
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larger scales for the dispersion function and their analogous features in the auto-correlation func- 
tion - are still present in Orion BN/KL (Figure [2]). However, already the lower resolution BEVIA 
observations, in particular at 1.3 mm, show larger irregularities than the high resolution W51 e2/e8 
data. Since they are probing different physical scales at the corresponding source distances (W51 
being 14 times further away than Orion), these larger irregularities might simply be due to different 
morphologies and structures in the magnetic fields. In all cases, the three lowest bins, showing a 
close to linear increase in dispersion, were again used to fit for the turbulent dispersion b. Values 
in the range of ~ 13° up to ~ 33° for the highest resolution are found (Tabled)). The corresponding 
turbulent to mean field strength ratios are around 0.2 for the BIMA and ~0.4 for the SMA obser- 
vations. Thus, over about an eight times increase in linear resolution, dispersion and turbulent field 
contribution increase by about a factor of three and two, respectively. 

All the auto-correlation functions have in common a rather sharp drop over the first three or 
four bins (Figure [2] right panels). With successively higher resolutions, this drop occurs within 
shorter scale ranges: within about 20" and 10" for the BIMA observations, and ~ 8" and less than 
5" for the SMA observations. In all cases the auto-correlation function regains after this first drop, 
either with a smooth slope (BIMA 3 mm, Figure [21 bottom right panel) or with several peaks. The 
characteristic large-scale polarization angle correlation length, equation ©, turns out to be in the 
range of ~ 66 mpc to ~ 15 mpc, which is a change of about a factor of four. 

The turbulent correlation function again shows a sharp drop over the first few bins. No sat- 
isfactory function was found for the BEVIA 1.3 mm observation. This is likely due to the incom- 
pleteness of the statistics (only few polarization segments) which the method in Appendix A relies 
upon. A, derived from the first few bins is between 36 and 9 mpc. 



5. Discussion 
5.1. Validity and Robustness of Approach 

Since the high resolution observations are revealing some differences and likely probe a dif- 



ferent magnetic field regime (section 15.41) . we re-assess here the validity of our approach when 
deriving the turbulent dispersion and turbulent to m ean field strength from e quation ©. In order 



to stay well within the linear range, the derivation in Hildebrand et al.l (|2009) assumes 6 < l k <C d. 
Whereas the turbulent magnetic field correlation length 5 is a fundamental limit related to the tur- 
bulence dissipation scale and the ambipolar diffusion scale L, the upper limit d, the typical length 
scale for variations in the large-scale mean field B , depends on the magnetic field under consid- 
eration. For the field threading a molecular cloud (envelope), d will be on the order of 100 mpc 
or more (e.g. envelope in W51 e2/e8). In the case of a collapsing core, the remaining large-scale 
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field is of the size of the core diameter, which is around 50 mpc for W51 e2 (|Tang et al.ll200 9b). 
It becomes clear that the allowable interval 5 < 4 <C shrinks with higher resolution. Neverthe- 
less, in such a cascading picture, a previously small scale field component becomes the large-scale 
mean field component at the next higher resolution. This still satisfies the assu mptions as long as 



the r esolution is not too high. It should be further stressed that the method in iHildebrand et al. 



(|20Q9f) is entirely independent of a mean field modeling, because of the restriction 4 *C d. Some 
values of the derived turbulent polarization angle correlation length A, (Appendix A) are compa- 
rable to the smallest measured scales or even larger. Therefore, X t ~ 8 < 4 is, strictly speaking, 
not valid everywhere. Should the correlation length indeed be of the order of 10 mpc, some of our 
observations are already resolving this scale. In such a case - referring to curve E in Figure 1 in 



Hildebrand et al.l (|2009l) - the turbulence contribution to the dispersion is probably underestimated 
because the fitting is done in a range where the turbulence contribution has not yet reached its 
maximum. The derived turbulent to mean magnetic field strength ratios are then lower limits. We 
note that the highest resolution SMA observation in Orion (Figure[2]top left p anel, 9 ~ 2.1 mpc) i s 



close to but no t yet re solving the ambipolar diffusion scale L (L ~ 1 mpc from lLi & Houd e (2008); 



Lazarian et all ((2004)). 



In order to probe the robustness of the approach, a threshold test is introduced. PAs with a 
continuum intensity abo ve a certain limit are excluded from the analysis. In the original work by 



Hildebrand et al.l (I2009Q only data with a magnetic field organization on large scales are analyzed. 
Our data set additionally contains data where the magnetic field is organized on smaller scales. This 
is for example the case for the hourglass-like pinched field lines in the co llapsing core of W 51 e2 



where the central pinched field lines correlate with the strongest emission (ITang et al. l2009bh . This 



threshold test addresses a possible concern whether there is a bias toward the strongest emission 
data. For Orion no relevant change in the dispersion function is found until discarding ~ 50% of 
the data (~ 70% cut in intensity). In particular, the turbulent dispersion derived from the fit shows 
a scatter within only ~ ±2° -3°, therefore not altering the original results and proofing them to be 
quite robust. This is verified for all the Orion data sets. Further excluding more than 50% of the 
data distorts the dispersion function, producing unreliable fits and turbulent dispersion values. The 
same results hold for W51 as observed with BIMA. The higher resolution observations of W51 e2 
and e8 however are found to be a little less robust. For both e2 and e8 a decrease of ~ -7° in the 
turbulent dispersion is found when excluding up to 30% of the data. Excluding more than 30% of 
the data leads to unreasonable results. In the case of a collapsing core (like e2 and possibly e8), 
gravity pulls in the magnetic field lines and presumably affects and dominates turbulence on some 
scales, bending the field lines on the shortest scales possibly more than what turbulence would 
do. Removing the strongest intensity data - which are in this case closest to the center of the 
collapsing core - is likely to reduce a possible gravity induced bias on the shortest scales. The 
observed decrease in the turbulent dispersion possibly reflects this. In the subsequent Figure [3j 
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these values are shown with a down-arrow (|) indicating this possible bias. 

Additionally, the dependence on the bin size has been checked. Decreasing the bin size by 
about 30%-50% of the synthesized beam leaves most of the results practically unchanged, with a 
typical ~ 1° to 2° down shift in the turbulent dispersion. This results from the denser sampling at 
the shortest scales which then extends the linear part of the dispersion function toward a slightly 
shorter scale. This in turn leads to a slightly lower value of the intercept at scale= 0. Whereas 
oversampling shows little effect, under- sampling with a twice as large synthesized beam signifi- 
cantly biases the analysis toward larger dispersion values. This can be understood from the steep 
slope over the first few bins. Averaging over large scales then increases the dispersion value. An 
exception to this conclusion are again W51 e2 and e8. Oversampling by a factor of two reduces 
the turbulent dispersion to ~ 15°. This is possibly due to the same reason as described above for 
the threshold test: Grouping the few values with a large dispersion in a separate bin will lead to a 
much lower dispersion for the remaining ones in another bin. Orion, from the BIMA 1 .3mm data, 
shows noticeable effects when the bin size is changed by only 10% to 20%. This is possibly due 
to the rather irregular dispersion function. 

In summary, the tests described above demonstrate that the analysis gives generally robust 
results with a small scatter. The possible exception is W51 e2 and e8, which might be biased 
toward too large dispersion values. In the remaining sections, the discussion will be based on the 
original values, with a reminder of the possible bias where necessary. 

The polarization angle correlation function < C(4) >, being mathematically related to the 
structure function via A0 2 (4) ~ 2(< C(0) > - < C(4) >), reflects the above discussed features 
analogously. We remark that, although this mathematical connection exists, the presented auto- 
correlation functions, < C(4) > and < Co(4) >, are calculated directly from the PAs which provides 
thus an independent consistency test. Mathematical relation and direct calculation are verified to 
lead indeed to the same results for < C(4) >• 

For the large-scale correlation function < Co(4) >, after separating the small scale turbulence 
contribution (appendix A), an immediate verification is not possible any more. No reference in the 
literature was found allowing a direct comparison here with other observations. Technically, the 
relative errors for the auto-correlation function are smaller than for the dispersion function (less 
than 1% compared to about 1 to a few percent), as a result of the error propagation of the sum of 
products compared to the sum of differences. 

Finally, we propose here to use the resulting large-scale polarization angle correlation length 
A as a quantitative measure for the typical scale d of variation in the ordered polarization structure. 
The derived values (Table [D match reasonably well with the above quoted empirical values which 
are estimated from the polarization maps. Our results are then based on the dispersion function 
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within 5 < 4 <C Ao. Having a statistical measure for d through the auto-correlation length A then 
also verifies the regime where the method based on the dispersion function (section 13.11 ) can be 
applied. Since A is an integrated and weighted measure, threshold test and bin size show only 
negligible changes. In the course of the tests described above, oversampling and under- sampling 
by a factor of two give variations in A of less than 10% for all observations. 



5.2. Comparison with Previous Results 



Generally, the dispersion functions derived from the present observation set for both W51 
e2/e8 and O rion BN/KL are less smooth than those derived from lower resolution single dish 
observations ( Hildebrand et al.ll2009 : Kirbyl2009 ). We speculate that these irregularities are tracing 
underlying changes and structures in the magnetic field morphology which become more manifest 
only at higher resolution (8), but are smoothed out at lower resolution. The previous larger scale 
(lower resolution) single dish observations show polarization patterns which extend to such radii 
where the magnet ic field straightens out as the gravitational influence is weak (e.g. OMC-1 in 
Schleuningl (11998b with the Kui per Airborne Observatory (KAO) at 100/iin with a resolution of 
~ 35" and in lHoude et aD hoo4) with the Hertz polarimeter at CSO with 6 ~20" at 350/im). M17 
SW, observed over a ~ 6' x 6' field with the KAO shows an overall ord erly field confi guration 
with a hint of being pulled into the cloud core by gravitational collapse (|Dotsonlll996r) . At the 
smallest separations, the average change in polarization is around 10°. T he sample of 12 Gal actic 
clouds observed with the KAO shows mostly organized large scale fields (|Dotson et al.ll200(X ). No 
dispersion functions were derived for this sample. 

Signs of a collapsing cloud in the center, with otherwise mostly straightened field lines, are 
observed for DR21 Main (IKirbyll2009l) . A smooth dispersion function is found with a dispersion 
of ~ 10° in the lowest bin and a linear increase to ~ 25° over the next three bins. The one case 
in which we observe an equally smooth dispersion is W51 observed with BIMA. Similarly to 
DR21 Main, this observation is tr acing the large-scale field envelope^ , but without yet revealing 
a collapsing core (ILai et al.ll200l|) . The lowest bin dispersion is again around 10°, then linearly 
increasing over the next four bins to about 30° (Figure Q] bottom left panel). The single dish 
(DR21 Main with CSO, spatial resolution ~ 0.2 pc) and interferometer (W51 with BIMA, spatial 
resolution ~ 0.1 pc) seem to be revealing comparable structures and dispersion values around 10° 
here. Such dispersion values then lead to a turbulent to mean field strength ratio of 0.1 to 0.15. 



The smaller scale field structures are likely to be beam-diluted because of averaging in the plane of sky, but not 
along the line of sight. In the opposite case, higher angular resolution data, e.g with the SMA, would not be able to 
reveal coherent structures on smaller scales. 
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These numbers are also consistent with the values reported recently in iHildebrand et al.l (I2009|) 
for OMC-1 and M17. Our remaining data are direc tly probing the coll apsing cores with clearly 
pinched magnetic field li nes for W5 1 (see Figure 6 in lTang et al.1 (|2009b() ) and a wrapped toroidal- 
like structure in O rion (ITang et al.l I2010T) . For this different physical regime - although previous 
observations exist ( Girart et al. 120061) - no analysis based on a dispersion function was performed 
which would allow a comparison. A turbulent to magnetic energy r atio of ~ 0.02 (turbulent to 
magnetic field strength ratio ~ 0. 14) is derived by Girart et al. d2006) based on a remaining mean 



dispersion after fitting a parabolic function to the field morphology. 

A further major difference between our data set (except W51 with BIMA) and the above cited 
previous works is the steeper slope (section 15 .31 Table Q]) over the first two or three bins. The 
dispersion increases typically by 20° to 30° or even more, compared to ~ 10° in the above cited 
cases. The values at the first bins are already around 20° or 30° whereas the other cases show 
numbers around 10° or less. The corresponding turbulent to mean field strength ratios are then in 
the range of ~ 0.2 to 0.54. 

Different methods, also aiming at constraining the turbulent to mean field ratio but not going 
through a structure function, have been explored by several authors. Based on an av erage PA 



disper sion S<p < 13°, a turbulent to magnetic field energy ratio has been estimated in lLai et al 



(120021) for NGC 2024 FIR 5 in the Orion B Giant Molecular cloud. Their energy ratio (<0.14) 
corresponds to a turbulent to mean field strength ratio <0.37. The ratio of mass-to-flux, M/$, 
is evaluated in iTroland & Crutcherl (|2008|) for a set of 34 dark cloud cores (~ 0.01 pc) from OH 
Zeeman observations with the Arecibo telescope. Their average ratio of turbulent to magnetic 
energy is ~2, the turbulent to mean field ratio therefore about 1 .4. Part of these data, in combination 
with GBT observations of the cloud envelope (~ 1 pc), were then used to p rovide support for a 
super- Alfvenic (weak magnetic field) turbulence model (ICrutcher et al]|2009l) . Average dispersion 
values (not dispersion fu nctions) around a mean magne tic field direction w ere also derived in 
Myers & Goodman! (II 99 1|) from optical polarimetry and in lNovak et al.l(|2009|) from submillimeter 
polarization. In the latter work the authors found a turbulent to mean field ratio in the range 
between 2.0 (intermediate field) and 0.52 (strong field) compatible with their observations. 



5.3. Comparison with Numerical Simulations 



The turbulent to mag netic field energy ratio is also investigated through numerical simula- 



tions. 



Ostriker et al.l (|2001l) analyzed the time evolution in a model cloud simulation with three 
different magnetic field strengths (strong, medium and weak field with (3 = c^/Bo/(4tt/5) = 0.01, 
0.1 and 1, respectively), but an identical initial turbulent velocity field. The perturbed mag- 
netic field energy reaches a maximum at about 0.1 - 0.2 times the Alfvenic crossing where it 
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accounts for about 20 - 50% of the total turbulent energy (kinetic energy and perturbed magnetic 
energy together). In their projected snapshot only the strong field model (B =14 /iG) leaves 
significantly correlated ordered polarization segments with a perturbed to mean magnetic energy 
fitwb - SB 2 /Bl « 0.27. This corresponds to a turbulent to mean field strength ratio of about 0.52. 
Although this is comparable to the numbers of our higher resolution observations around collaps- 
ing cores, it is not obvious to match their snapshot in time evolution within our sequence of low 
and high resolution data. 

On the other hand, recent simulation results with supersonic and super- Alfvenic turbulence, 



in combination with simu 



ated Z eeman measurements, find (3 tur b ^ 1 for cores with a radius 



0.2-0.8 pc (ILunttila et al.l 120091). Th i s seem s to be in favor of the turbulent to magnetic energy 
ratios observed in lTroland & Crutchen (I2008|) . who found an average /3 turb w 2 for comparable core 



sizes and source distances. These findings support a star formation theory with super-Alfvenic 
turbulence. 



Related to f3 tur b, PA structure functions (of second order), A0 2 , are analyzed in lFalceta-Goncalves et al 



(I2008|) with the goal of discriminating between sub/supersonic and sub/super- Alfvenic models. In- 
dependent of the magnetic field orientation with respect to the line of sight, sub-Alfvenic models 
tend to have a power law index a ~ 0.5, A0 2 (4) ~ 1%, whereas super-Alfvenic models show a 
flatter slope with a ~ 0.3. For the set of our observations, a (derived from the first two bins for 
W51 and three bins for Orion, Tabled)) is in the range of 0.3 to 3.2. Although this is closer to 
sub-Alfvenic models, the rather steep slopes in our structure f unctions followed by a pl a teau m ake 
a direct comparison with the smoother structure functions in iFalceta-Goncalves et al.l (J2008) not 
obvious. 

In summary, numerical simulations provide results for the turbulent to magnetic field energy 
ratio and the slope of the polarization structure function for a series of different models. Some 
are in agreement with our findings. For a detailed comparison, the difficulty lies in matching 
observational resolution and star formation stage with the model time evolution in the simulation. 



5.4. Dependence on Physical Scale 

The angular resolutions obtained from the BIMA and SMA observations vary by about a fac- 
tor of ten, spanning a range in physical scales at the source distances from ~ 70 mpc to ~ 2. 1 mpc 
(TableQ}. As discussed in sectionH the resulting dispersion functions (and auto-correlation func- 
tions) show some common tendencies but differ in their detailed characteristics and numbers. Here, 
we address the question whether the subsequently higher resolution observations reveal a depen- 
dence on physical scale. This question is also motivated by the observed polarization maps where 
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an envelope on the largest scales (BIMA. lLai et al.l (1200 1|) ) and two collapsing cores wit h a possi 



ble hou rglass -like magnetic field morphology on the smallest scales (SMA, Figure 6 in lTang et al 



(I2009b|) ) are found in the case of W51 e2/e8. For Orion, the lowest resolution BIMA data and the 
highest resolution SMA data again reveal very different magnetic field structures. Consequently, 
we argue that these polarization observations are indeed probing different regimes in the star for- 
mation process. The fact that increasingly higher resolutions still reveal detectable local coherent 
structures, very likely means that large-scale structures are smooth enough for the local structures 
to be distinguishable. This must hold true for both along and across the line of sight, because 
otherwise features would be washed out. 

Figure [31 top panel, shows the tendencies for the turbulent to mean field strength ratio as a 
function of the physical scales, as derived from the Figures Q] and [21 As relevant physical scale the 
synthesized beam re solution is assumed. In addition to the BIMA and SMA data set, the lower 
resolution data from Hildebrand et al.l ((2009) for Orion, DR21 and M17 are also added. Over the 
observed scale range (~ 70 mpc to ~ 2. 1 mpc), the turbulent to mean field strength ratio increases 
with smaller scale by about a factor of ten. 

Further investigating this apparent dependen ce on phys i cal sc ale in the top panel in Figure [31 
w e discuss a possib l e be am resolution effect. iHoude et al.l d2009|) . in an expansion on the work 



in 



Hildebrand et al.l (|2009|) . considered the signal integration aspect in their analysis. They find 
that the turbulent component of the dispersion function at scale zero, b 2 (0), is then the square of 
the turbulent to large-scale magnetic field strength divided by the number of independent turbulent 
cells N probed by the observation: b 2 (0) =< Bj > jB\ ■ l/N. As further derived, N is directly 

/ jC"2 W7 2 \ /\ / 

a function of the beam size and can be written as N = —, — 4 — , where 5 is the turbulent field 
correlation length, A' is the effective depth of the molecular clo ud along the line of sight and W 
is the beam radius. The equation for /V in IHoude et al.1 d2009|) is derived assuming a (circular) 
Gaussian beam and Gaussian turbulent auto-correlation functions. For a given source, the smaller 
the beam size, the smaller /V will be, and the bigger b will be if the ratio < Bj > l l 2 /B remains 
constant, or if at least 1 / \/N grows faster than < Bj >^ 2 /B decreases. Different beam sizes can 
therefore mimic a trend in the turbulent to mean field ratio. The top panel in Figure [3] apparently 
confirms this expectation. In the following we aim at revising the top panel in Figure [3] by taking 
into account the beam resolution effect, in order to reveal the ne t change in t he tur bulence to 
mean field strength ratio over scale. Starting from equation (44) in IHoude et al.l (|2009|) . the scale 
independent term describes the intercept at scale zero in their Figure 2. This is equivalent to 
our dispersion at scale zero, but with an additional factor containing the beam integration effect. 
We, therefore, rewrite the turbulent to mean magnetic field ratio as: < Bj >^ 2 /B = a//V ' / 'v2- < 
A(p 2 >|io ' where < A0 2 >){^ is the turbulent dispersion derived from our fitting results in Figured] 
and|2] A' is approximated with the size of the most extended detected structure, whi ch is roughly 
the size of the maps in the SMA and BIMA observations. <5(~ 16 mpc) is adopted from lHoude et al 
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(|2009|) for Orion. The beam radius W is derived from the synthesized beams of our SMA and 
BEVIA observations. The estimated number of turbulent cells is listed in Tabled] Figure [31 middle 
panel, shows the resulting beam corrected turbulent to mean magnetic field strength ratios. After 
correction, both Orion and W51 show a close to constant ratio over scale, with a slight trend of 
a decreasing ratio toward smaller scales. Whereas the ratios for Orion are in the range between 
0.30 and 0.44, the ratios for W51 (0.7 to 1.27) are around the equipartition limit of turbulent and 
magnetic field strength. This might be an indication that the role of turbulence and magnetic field 
changes over scale. In the case of W51, there is possibly a transition from a turbulence dominated 
to a magnetic field dominated scenario. 

Taking into account the beam integration modifies the picture first presented in the top panel. 
Differences in measured values of < A</> 2 >^ can result from a combination of a beam integration 
effect and a variable turbulent to mean field ratio. Both contributions can vary in their significance 
with scale. Correcting for the beam integration effect reveals the net change over scale in the turbu- 
lent to mean field ratio. The beam correction (~ y/N) is of importance for lower resolution (larger 
beam) observations, but becomes less important for higher resolution observations approaching 
the turbulent cell size. Intuitively, one might expect a turbulent to mean field strength ratio which 
decreases with smaller scale as the turbulence dissipation scale is approached. The observed trend 
seems to be in support of this, revealing the beginning of a slight decrease in the ratio. An even 
higher resolution observation at the mpc scale or smaller might then reveal a breakpoint (turn over) 
in the scaling. On the other hand, the panels in Figure[3]only display a ratio. It is still possible that 
the turbulent field strength increases with smaller scale, but at a slower rate than the mean field 
strength. A power law, oc P, is fit to the Orion and W51 data. Table [2] summarizes the results for 
both the uncorrected and the beam corrected scalings. 

We finally remark that there is some uncertainty left in the beam correction due to our approx- 
imation for A'. Should the detected cores have structures along the line of sight which significantly 
differ from the detected extension across the plane of sky, the number of turbulent cells N could 
vary. Underestimating /V by a factor of two would change the ratio by y/2. In order to align Orion 
and W5 1 , N would need to be underestimated or overestimated, respectively, by about a factor of 
four. Similarly, a correlation length 5 twice as large as in Orion would reduce the ratios for W5 1 
to approximately the ones found for Orion. However, remaining uncertainties in both A' and 5 do 
not change the slope of the scaling, unless they significantly vary with resolution. Thus, within 
these uncertainties Orion and W5 1 both show indications of a decrease in the turbulent to mean 
field strength ratio over scale, but with a difference in scale specific ratios possibly reflecting the 
source environment. 



The scaling of the large-scale polarization angle correlation length A shows a close to straight 
line (Figure |3j bottom panel), with only a small difference between W51 and Orion. A decreases 
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by about a factor of fifteen over the sampled physical scales. We remark that the presented values 
for A are possibly lower limits, because of the relatively small fields of view sampled in our 
observations. Since A is a normalized measure, equation ©, beam integration effects are likely 
to cancel out, unless they additionally depend on the spatial scale l k . The turbulent correlation 
length X t is presumably independent of scale. This is not fully verified from our data set (Figure 
|3l bottom panel). In some cases this is due to incomplete turbulence statistics, which then do not 
allow to set a reliable cut off criteria. Additionally, a certain beam integration effect is probably 
left in V(A(j) k=l ), equation ®, which leads to broader or narrower turbulence distributions. This 
again affects the cut off criteria through //a (Ap pendix A). For com parison, a turbulent magnetic 
field correlation length of « 16 mpc is derived in lHoude et al.1 (I2009Q for OMC-1. 



6. Summary and Conclusion 



A set of interferometric polarization observations with resolutions in the range between 1" 
and 0."7 (~70 mpc to ~ 2.1 mpc) is analyzed to derive statistical properties of magnetic field 
and turbulence from large to small scales during the star formation process. Our data set covers 
structures from the large-scale cloud envelope to the collapsing cores. The highest resolution data 
are close to the expec ted ambipolar diffusion scale (~ 1 mpc). We apply and expand the method 
developed recently in lHildebrand et al.l (|2009|) . The main results are: 



1. The turbulent field dispersion shows a steeper slope and larger values at the shortest scales 
with increasing resolution. Accordingly, the resulting turbulent to mean magnetic field 
strength ratio increases with smaller scale over the entire range in physical resolution. This 
is without taking into account a beam integration aspect, and is in agreement with earlier 
theoretical expectations. 

2. The sequence of low and high resolution observations does not only zoom in onto the same 
magnetic field structure, but it is probing different morphologies and different stages in the 
star formation process. This is also supported by the polarization maps. When taking into 
account a beam integration aspect, both Orion and W51 show a close to constant turbulent to 
mean field strength ratio over scale, with a slight trend of a decreasing ratio toward smaller 
scales. Whereas the ratios for Orion are in the range between 0.30 and 0.44, the ratios for 
W51 (0.7 to 1.27) are around the equipartition limit of turbulent and magnetic field strength. 
This might be an indication that the role of turbulence and magnetic field changes over scale. 
In the case of W51, there is possibly a transition from a turbulence dominated to a magnetic 
field dominated scenario. Our observation set therefore also provides information for the 
time and spatial evolution of these quantities during the star formation process. 



-18- 



3. Based on the polarization angle correlation function a characteristic large-scale correlation 
length A is defined. This can be used as a quantitative criterion to define the scale over which 
the mean polarization structure varies. A decreases with higher resolution. Additionally, 
starting from an ensemble of measured polarization position angles, a method is proposed to 
separate statistically large-scale from turbulent contributions. This leads to a definition of a 
turbulent correlation length \. 

The authors wish to thank the referees, Roger H. Hildebrand, Martin Houde and John E. 
Vaillancourt for their comments and explanations which provided important further insight. 

A. Polarization Angle Correlation Function 

An observed position angle PA t = 0,- is the superposition of a large-scale polarization structure 
and a smaller scale turbulent component. Consequently, a correlation length directly derived from 
a measured ensemble of 0,- contains both small scale and large-scale correlations mixed. In order to 
calculate a large-scale (mean) polarization angle correlation length separately, one needs to isolate 
the mean contribution 0; O and the turbulent contribution dept. In general, this will not be possible 
for a single or only a few PAs, unless one makes very specific assumptions as, e.g. a given uniform 
magnetic field direction. However, it is, as outlined in the following, possible to separate the two 
contributions in a statistical way. 

Splitting each PA t into a large-scale and turbulent part, ; = 0, o + Sfa, the correlation product 
for a scale k is written as 

^2 & ' ~ S ' + 2 S ^° ' + ^ ' ^ ( 6 ) 

where the mixed correlations 0,- o • S<f>j and 7 - o ■ 50; are set identical when evaluated within the 
same scale range. The summation extends over l k < r t j < 4 +1 where = sj (jc,- - xj) 1 + (y,- - yj) 2 
with coordinates (x,,y,) for the position angle </>,-. The normalization factors are omitted here, and 
reintroduced later (therefore the sign '~'). Since observationally only 0, is accessible and we wish 
to isolate 4>io ' 4>jo, we replace 0,o = 4>i — 54>i in the mixed correlation, which leads to: 

^0;o-0yo~ ^2<pi-(pj-2^2 4>i -tyj + ^tyi ■ 5(f) j. (7) 

In order to further proceed, the turbulent contribution 5(f) needs to be quantified. Although 
84>i is not directly observable for an individual position angle 0,, its distribution V(5(p) can be con- 
structed. At the smallest observable scale k = 1 , the difference in two position angles, A0, ; = 0, — 0y, 
is mostly reflecting the difference in their turbulent contributions. This is similar to the turbulent 
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dispersion value whic h is derived by ext r apolat ing from the smallest scales to the intercept at scale 
k = (section 13.11 and iHildebrand et al.1 ([2009J)), with the difference that Afaj can be positive or 



negative. The turbulent distribution V(5fa) can then be constructed as 

W) = \ b -^m • ^( A ^=i) (8) 

where the distribution of Afa = \ is evaluated from the derived differences in position angles at the 
smallest scale k=l, and the factor 1/2 results from assigning half of the difference as turbulent 
contribution to each of the two PAs. The additional factor bj < Afa 2 >^ down- weights the mea- 
sured differences at scale k = 1 to the limiting scale k = 0, based on the result from the fit for the 
turbulent dispersion function (section I3TI) . 

Typically, the turbulence distribution V(5fa) is expected to be Gaussian around zero. This 
further means that for a large enough, statistically complete sample, the second term on the right 
hand side of equation © will vanish due to the symmetry of T>{5fa). In order to verify this, the 
mixed term, fa ' 1S calculated by randomly choosing a 5(f) j value from T>{5fa). Except for 
the smallest scale, there are fewer correlation products than sample values 5(f)j in T>(8fa), and the 
mixed term will therefore typically not converge to zero. This statistical sampling limitation can be 
overcome by repeating the calculation and averaging the results. Typically, for 100 runs or more the 
results converge and then really probe the statistical completeness of the turbulence distribution. 
Figure |4] illustrates this for the case of W51 observed with BEV1A: V(5fa), top left panel, is very 
close to a Gaussian distribution around zero (fi ~ 0.9°). The resulting mixed term ^2 fa ■ 5(f) j, the 
separation between the dotted and dashed line in the top right panel, is close to zero. 

It remains to evaluate the third term on the right hand side of equation ©, Y 5fa ■ 5(f>j, which is 
the small scale turbulent correlation. Whereas the mixed term extends over all scales, because it is 
a correction term due to a local turbulent fluctuation with any possible position angle in each scale 
range, the turbulent correlation is exp ected to extend ove r a spatially limited area, characterized 



by the turbulent correlation length. In iHoude et al.l (|2009h this is taken into account by assuming 



a Gaussian window function where the width is the turbulent correlation length. In their further 
analysis the correlation length is then a fitting parameter. Here, we are directly making use of 
the turbulence statistical distribution to set a cut off in scale for the correlation term ^,6 fa ■ S(j)j. 
As in the case of Y fa ■ (f>j, the relevant quantity for the correlation is the difference, A = \Sfa — 
S(f)j\, between a pair of turbulence fluctuations 5 fa and 8(f) j. The distribution of A can again be 
constructed by calculating differences between two randomly chosen values from the distribution 
V(5fa). V(A) is shown for W51 (BIMA) in Figure HI middle left panel, with a resulting mean 
difference in turbulence fluctuations //^ ~ 3.8°. We are adopting as a cut off measure. This 
means that 5 fa ■ 5(f) j is counted as a correlation pair for a a certain scale k if the difference in the 
corresponding fa ■ (f>j pair, A0 = \ fa — (f»j\, is smaller than [1a- Otherwise the correlation in fa ■ <f>j 
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is considered to be dominated by the large-scale magnetic field. Alternatively, the cut off criteria 
can be refined by calculating the probability function of A0 being of turbulent origin: to that 
purpose, the normalized distributions of P(A) and Z>(A0) are weighted against each other, Figure 
HI bottom left panel. (When working with the normalized distributions, the weighting factors (see 
below) are not yet taken into account.) Comparing the distribution V(A) with the distribution 
V(A(p) provides thus a tool to separate small and large-scale contributions. Both the /iA cut off 
and the more sophisticated weighting with the probability function lead to very similar results. We 
therefore adopt the simpler criteria with //a m the following. Although T>(A) is only known in 
a statistical way, with the same /^a used for each correlation pair, this average value can still be 
used to check against the measured difference in each correlation pair 0, • (fij. Technically, we thus 
evaluate: 

6<f>rS<f>j = <f>r<f>j if A0 < /i A (9) 

When adding the two correlation parts in equation C7]), Yl4>i' 4>j an d ' a P ro P er nor- 

malization needs to be introduced. A factor l/(Nk+Ng,k) with N$,k counting the correlation pairs 
satisfying equation ©, normalizes the sum to one. The large-scale correlation function, < Co >, 
for each scale k can then be written as 

with the condition in equation ©. 

The isolated turbulent correlation, 1 /N$ t k • ^2 $4>i • 84>j> ls shown in FigureHl middle right panel. 
A weighting factor A^/A^ is taken into account to directly compare its statistical significance with 
< C > k in the right upper most panel (weighting factor A^/A^ = 1). Thi s additional weighting 



provides a spatial filter similar to assuming a Gaussian window function (|Houde et al.ll2009|) . It 
acts like a probability, for each scale k, that a measured A</> < /iA contributes to the turbulent 
correlation. We note that this is also valid when adopting the probability function as a cut off 
criteria, as mentioned above. 

The bottom right panel in Figure 0] shows the final large-scale correlation function < Co > k 
by combining the upper two panels with the proper normalization from equation (flOl) . Although 
the turbulent correlation shows a relatively sharp decrease over the smallest scales, its overall 
correction to < C > k (dashed line) is small. This is due to the small statistical significance (in terms 
of number counts) of the isolated turbulence differences in PAs. Adding up the various correlation 
factors with their statistical weight is one possible choice. Its most immediate advantage is that it 
gives directly the correct normalization. 

Finally, we note that the characteristic large-scale polarization angle correlation length, Ao, is 
thus virtually unchanged compared to when it is directly calculated from the measured ensemble 
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of 0,-. From the middle panel, additionally, a turbulent polarization angle correlation length, X t , can 
be estimated in an analog manner. 
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Fig. 1. — Left Panels: Dispersion (square root of second order structure function), binned to multiple integers of 
the resolution of each observation for W51 e2/e8. The assigned scale for each bin is obtained from averaging all 
the scales within a bin, which is very close to the synthesized beam in the first bin and close to the center for all 
the following bins. Correlated data points below the synthesized beam resolution are removed. Error bars are very 
small at the smallest scales due to the sample variance factor (large number of data points), and grow for the larger 
scales. The red line is the fitting result including the first two or three bins following equation where the turbulent 
dispersion b is found from the interception of the fit at scale=0. Right Panels: large-scale polarization angle correlation 
function <Cq> (solid line) binned to multiple integers of the resolution of each observation. The weighted integral of 
the curve measures the characteristic large-scale polarization angle correlation length Ao. The number of data points 
(pairs of PAs) within each bin is displayed (□) with the axis on the right hand side of each panel. The same number 
of data points are also used for calculating the structure function, except where correlated data points are removed. 
Binned values are connected with straight segments for visual guidance only. The dashed line shows the turbulent 
polarization angle correlation function, separated by the method described in Appendix A. The dotted line is the raw 
auto-correlation function < C >, directly derived from a measured ensemble of PAs without separating large and small 
scales. At the distance of W51 (~ 7 kpc), 1" corresponds to ~ 30 mpc ss 6190 AU. 
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Fig. 2. — Identical to Figured! but for Orion BN/KL. No satisfactory solution was found for Orion 
(BIMA 1.3mm). The correlation shown is the raw polarization angle correlation function < C >. 
At the distance of Orion (~ 480 pc), 1" corresponds to ~ 2.3 mpc 470 AU. 



NO BEAM INTEGRATION CORRECTION 





Fig. 3. — Top panel: turbulent to mean field strength ratio < Bj > 1//2 /Bo from various observations as a function 
of physical scale, as derived from the Figure [TJ and [2] As relevant physical scale the achieved resolution (synthesized 
beam) a t the source distance is a ssigned. Added to the BIMA and SMA data set are DR21 (*), M17 (□) and Orion 
(+) from Hild ebrand et al. (2009). The down arrows for W5 1 e8 and e2 mark the values after correcting for a possible 
gravity induced bias (~ -7°) as discussed in section IBTT1 The middle panel shows the same ratio, but taking into 
account a beam integration correction as described in section 15.41 Resulting power law fits (oc P) in both panels 
are shown separately for the original W5 1 data (blue, dashed line), Orion (red, dot-dashed line) and the entire data 
set (black solid line, only top and bottom panel). The indexes are summariz e d in T able |2] For illustration shown 
is also the expected ambipolar diffusion scale at ~ 1 mpc from iLi & Houdd J2008I) (dashed vertical line) and the 
turbulent - magnetic field strength equipartition line (dot-dashed line). Bottom panel: large-scale polarization angle 
correlation length, Ao, as a function of physical scale, together with the turbulent polarization angle correlation lengths 
("diamonds^. 
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Fig. 4. — Illustration of the analysis of the polarization angle correlation function for W51, BIMA. Top left panel: 
the turbulence distribution T>(Scf>), equation ©, with mean n « 0.9. The resulting correction to the raw auto-correlation 

< C > - ideally identical to zero for a Gaussian turbulence distribution around zero - is shown in the top right panel 
with the dotted line. The dashed line shows < C > calculated directly from the measured position angles r/>,. Middle 
left panel: the distribution of differences between two turbulence values, 2?(A), derived by random sampling two 
values from the distribution in the upper left panel. The larger number counts results from oversampling (calculating 
more difference pairs than combinations in T>(Scf>)). This also ensures a smooth distribution with a mean converging 
to fi a ~ 3.8 (dashed line), which is used as a cut off criteria for the turbulent polarization angle correlation function 
(middle right panel). A statistical weight (number counts) is applied (see appendix). The large-scale auto-correlation 

< Co > (bottom right panel, solid line) is obtained by adding < C > and the turbulent correlation with the proper 
normalization from equation (11 Ot . The dashed line shows again < C > for comparison. The only small correction 
from the turbulent correlation results from its small statistical weight compared to the full distribution of differences 
in PAs which is used for < C >. Bottom left panel: the normalized distributions of 2?(A) (blue line) and 2?(A0) (red 
line) which are used to derive the probability function (black dashed line) as an alternative to the cut off criteria in the 
middle left panel. 
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Table 1 . Observation Summary and Results from the Statistical Analysis 



w51 



Observation Analysis 



<s 2 > l/2 / <B 2 > 1/2 ' 

A d ref b — ^ — m a A A, N ( — ^ — 

(mm) (") (mpc) (mpc) (deg) (deg/") (mpc) (mpc) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 



BIMA 


1.3 


2.3 


69 


300 


SMA - e2 


0.87 


0.7 


21 


60 


SMA - e8 


0.87 


0.7 


21 


60 



Orion 



BIMA 


3 


7.0 


16.1 


69 


BIMA 


1.3 


3.4 


7.8 


69 


SMA" 


0.87 


2.8 


6.4 


23 


SMA'' 


0.87 


0.9 


2.1 


23 



I 


6.4 


0.08 


2 


II 


27.6 C 


0.36 


35 


II 


54.0 C 


0.89 


10 


III 


14.6 


0.18 


3 


III 


12.9 


0.16 


10 


IV 


24.7 


0.32 


7 


IV 


33.1 


0.45 


10 



1.1 


230 


45 


122* 


1.22 


1.3 


73 


25 


4 d 


0.70 


0.3 


63 


30 


4 d 


1.27 


1.6 


66 


36 


6 


0.44 


3.2 


46 




4 


0.30 


1.3 


35 


10 


1 


0.33 


0.7 


15 


9 


1 


0.43 



Note. — All statistical quantities are obtained from the Figures [T]and|2] with the method described in section]!] 
'''observing wavelength 

' 2 ' synthesized beam resolution, also used as the separation in between the bins \; for elliptical beams the geometrical mean is adopted, except for Orion 
BIMA 3mm where the semi-major axis is used because of its very elliptical beam 

'''synthesized beam in physical scale at source distance 

' 4 ' scale of variation in the large-scale polarization structure, empirically estimated from polarization maps by visual inspection 
^'references: (I) lLai et alJkoOlh . (II) lTang et all l2009bl) . (III) lRao et all Il998h . (IV) lTang etallfcOld) 
' 6 'turbulent magnetic field dispersion 
' 7 'turbulent to mean field strength ratio 

' 8 ' slope over the first two bins (W51) or three bins (Orion) in the dispersion function < A</> 2 > 1/2 

'''power law index of the second order structure function, < A0 2 >~ l a , over the first two bins (W51) or three bins (Orion) 
''"'plane of sky projected large-scale polarization angle correlation length 

'"'plane of sky projected turbulent polarization angle correlation length, derived from the method in Appendix A 
<12 'number of turbulent cells contained within telescope beam 
' 13 'turbulent to mean field ratio corrected for beam integration effect 
"compact and subcompact observation combined 
b compact and extended observation combined 

c possibly influenced by gravity, with a bias of ~ —7 deg (section l5~Tl 

d based on approximating A' (effective depth of the molecular cloud along the line of sigh t) with the size of emission which is roughly the size of the 
maps in the SMA and BIMA observations. The turbulent correlation length 5 is adopted from lHoude et alj ilQO^D . 
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Table 2. Power Law Indices 7 for Scaling Relations 



parameter 


W51 


Orion" 


combined 6 


< B] > 1/2 /B Q (uncorrected) 


-1.51 


-0.43 


-0.33 


(< B 2 t > 1/2 /B ) N (beam corrected) 


0.17 


0.07 




Ao 


1.02 


0.76 


0.75 



Note. — All power law indices 7 are derived by fitting the relation 
oc P as a function of physical scale / (Figure [3]). Separate fits are 
performed to subsets (W51, Orion) and the entire data set. 



includes Orion from iHildebrand et all (120091) for < B] >'/ 2 /B 
(uncorrected and beam corrected) 

b W51 and Orion c ombined. Orion, Ml 7 and DR21 from 
Hildebrand etaP too4) are included for < Bj >^ 2 /B (uncorrected). 
No joint fit is attempted with the beam corrected data. Ao is fit from 
the combined SMA and BIMA data, analyzed in the work here. 



